Unveiling the dynamics of the universe 
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We present two methods for describing, from a pedagogical point of view, the solutions of Ein- 
stein's equations applied to a homogeneous and isotropic universe. In the first method, we define an 
effective gravitational potential of the universe, and the second method makes use of the fact that 
we can define a dynamical system of equations. The methods are applied to different cosmological 
models whose properties are discussed in turn. The ultimate intention is to provide simple examples 
to be revised in a first Cosmology course for undergraduates. 

PACS numbers: 98.80.-k,98.80.Jk 



I. INTRODUCTION 

In any introductory course on Cosmology, some ques- 
tions undoubtedly arise once the students become famil- 
iar with the equations describing the evolution of a ho- 
mogeneous and isotropic universe. These are: a) which 
is the universe made of?; b) what is the curvature of the 
universe? c) how did the universe expand in the past 
according to the answers in a) and b)?. 

It is interesting to note that present technology has al- 
lowed the humankind to partially answer those questions 
with some important accuracy[l|,|^|3]- However, such an 
improvement in our knowledge of the universe will take 
still some years to reach the textbooks for undergradu- 
ates. 

In this paper, we briefly review the basic equations 
of the evolution of the universe, and describe how the 
latter evolves according to its material contents and cur- 
vature. It is not our aim to give an exhaustive study of 
the gravitational dynamics that arises from the equations 
of Einstein's General Relativity (GR), but rather to give 
some tools to help students to understand the richness 
of the solutions of the Einstein's equations in the case of 
a homogeneous and isotropic universe. In other words, 
we will focus our attention in the answer to question c) 
above. 

In the same line, our intention is also to provide simple 
pedagogical exercises, and to solve them with the help of 
analytical methods which are of widespread use in the 
specialized literature. The chosen examples in this paper 
are appropriate, in our opinion, for a first undergraduate 
course on Cosmology. 

We shall not discuss how we can determine the mate- 
rial content of the universe and its true nature through 
observations, nor how its curvature is measured. We will 
only mention briefly how such quantities are obtained, 
and give the reader some interesting references and inter- 
net links where more detailed information can be found. 
However, non-experts will find a comprehensible sum- 
mary of modern Cosmology inQ, and some interesting 
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questions reviewed and answered in0,|^. 

The present manuscript is organized as follows. In 
Sec. m we review the metric quantities that describe 
a homogeneous and isotropic spacetime, and how they 
are influenced by the material content of the universe 
through the equations of Einstein's General Relativity. 
We describe the types of matter we shall consider, and 
how each one is classified according to their equation of 
state. 

We shall introduce the so-called density parameters, 
which measure the relative contribution of each compo- 
nent to the total material content of the universe. These 
parameters will play a central role in the subsequent cal- 
culations. Due to its particular importance, we will take 
some space to discuss the significance of the curvature's 
density parameter. 

In Sec. mil we describe how the Friedmann equation 
can be used to define an effective gravitational potential 
of the universe, so that we can visualize the expansion 
of the universe in a similar fashion as one describes the 
motion of a single particle in Classical Mechanics. 

In Sec. IIVI we use again the Friedmann equation and 
the equations of motion of each single fluid to obtain a 
set of differential equations that allows us to see the ex- 
pansion of the universe as the solutions of a dynamical 
autonomous system. The latter formalism will be used 
to see whether there is any attractor behavior in the cos- 
mological solutions. 

The formalisms developed in the previous two sections 
is applied to three particular examples in Sec. These 
are: the actual standard cosmological model, also known 
as the Concordance Model (CM); Einstein's static model 
of the universe; and a CM with an arbitrary component 
of what we shall call dark energy. 

Finally, Sec. IVII is devoted to conclusions and general 
comments. 



II. MATHEMATICAL BACKGROUND 

We start with the so-called Cosmological Principle 
(CP), first proposed by Einstein, which states that the 
universe we live in is homogeneous and isotropic in large 
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scales. By large scales we mean scales much larger than 
the size of a typical galaxy. 

The CP is a working hypothesis, and it is the sim- 
plest assumption we can make about the properties of 
the spacetime of the universe as a whole. The interested 
reader can find a more detailed discussion on the histor- 
ical and philosophical relevance of Einstein's CP in|^,Q. 



A. Space geometry 

The remarkable thing is that the CP suffices to fix 
the metric of the spacetime a homogeneous and isotropic 
universe must have. It is called the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric, a metric with con- 
stant curvature, whose line element is usually written as 
(in units with c — 1) 



-dt' 



dr^ 



1 — kr"^ 



(1) 

where the coordinate t is called the cosmic time, a(t) is 
the (time-dependent) scale factor, and k is the curvature 
constant. 

The spatial part of the metric is written in terms of co- 
moving coordinates, as they remain fixed as the universe 
expands. Actually, the physical and comoving distances 
at a given time t are related through Xph = a{t)xcom- 

Notice that dO. is the usual spherical solid angle ele- 
ment, and r is a comoving radial coordinate. In general, 
for a fixed r and t, we see that the circumference of a 
circle is equal to C = 27rr, and the surface of a sphere is 
S — 47rr^j3l. But, the proper radius of these objects, for 
fixed 9 and (p, is given by 



(.= 



dr 



■3= sin ^{^/kr 



r 

■7= sinh" 



for A: > , 
for A: = , 
'\^/^r) for A: < 0. 

(2) 

As a consequence, we see that for for fc = 0, we recover 
the usual Euclidean relations C/£ = 27r and S/i'^ — An. 
However, if the curvature constant is positive fc > 0, then 
C/£ < 2tt and S/£-^ < An. The opposite happens for 
a negative curvature fc < 0, and hence C/£ > 2tt and 
S/i"^ > An. 

The expansion of the universe is encoded in the time 
derivative (which we denote by a dot) of the scale factor, 
with which one defines the Hubble parameter H = a /a, 
named after the astronomer Edwin Hubble. Its impor- 
tance comes from the fact that nearby galaxies recede 
with a velocity proportional to their (physical) distance 
from us, as ?; = ffccph [gllTollTTf . 



B. Cosmological equations 

The equations of motion of the universe are given by 
Einstein's GR, which is our current fundamental theory 



of gravitation|l[7l li. ITlllllTlIT^ . For this, we need to 
specify the kind of matter the universe is made of, and 
give its corresponding energy-momentum tensor. The 
usual choice is that of a homogeneous and isotropic per- 
fect fluid, which is well described and defined only by its 
energy density p{t) and its isotropic pressure p{t). 

There are two equations coming from GR, relating the 
geometry and material contents of the universe, and there 
is another equation comiiig from the conservation of the 
energy-momentum tensor 0, 0, 0, 0| . For a heuristic 
approximations to the cosmological equations see[l5l[l6|. 
The equations of motion are 



H = -AnG^{p^+p^ 

i 

Pi = -3H{pi +p^) , 



(3a) 
(3b) 



together with the so called Friedmann (constraint) equa- 
tion 



8n G fc 



(4) 



We have allowed for the existence in the universe of 
more than one perfect fluid; with the explicit restriction 
that they do not interact with each other except gravita- 
tionally. This means that each of the perfect fluids obeys 
a separate conservation equation, namely Eq. Ij3b|l . 



C. Types of matter 

Only (n + 1) of the above {n + 2) equations arc in- 
dependent, and we need (n) extra equations in order to 
solve for the {2n+l) unknowns a{t), pi{t) andpi(t). The 
needed equation is what is called an equation of state 
relating the energy density and pressure of each perfect 
fluid, we will follow the common wisdom and assume a 
barotropic equation of state in the form 



Pt = hi - l)p» : 



(5) 



where 7, is the equation of state of the i-th fluid. In this 
case, Eq. (|3b|) integrates to 



(6) 



where p^.o is an integration constant. The most usual 
values of 7^ are: 

• Radiation and relativistic particles, 7^ — 4/3 and 
Pr ^ a ^; 

• Pressureless matter (dust or non-relativistic parti- 
cles), 7™ = 1 and p„i a~^; 

• Cosmological constant (A or vacuum energy) , 7a = 
and PA ~ const 
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y However, it is widely accepted that the allowed range 
can be 2 > 7 > 35j. 

Eqs. and Q shows that the curvature of the uni- 
verse contributes to the dynamics. Actually, we can 
think of curvature as a kind of special perfect fluid with 
7fc — 2/3, so that pk ^ a~^, as it is indeed the case. 



D. Critical energy density and density parameters 

The Friedmann equation can be rewritten in the 
dimensionless form 



^ Pc Pc ^ 



(7) 



We have defined in here the critical energy density pc and 
the density parameters fli as follows 



8nG 



n, = 



-k 



(8) 



The density parameter for the curvature component flk 
is defined together with the negative sign, so that it is on 
equal footing with respect to the density parameters of 
the material components l^j . 

It can be seen from Eq. that if the whole material 
content of the universe equals the critical energy density, 
= Pc then the curvature of the universe should 
be zero; that is, it should correspond to a flat universe. 
Likewise, if J^i Pi > Pc Pi < Pc) then the universe 
has a closed (open) spatial geometry. 

The critical energy density depends on the value of the 
Hubble parameter, and thus is a time-dependent quan- 
tity. However, the current value of pc,o can be measured 
directly from the redshift of nearby galaxies and other 
objects, seeP, for recent results. 

On the other hand, each of the density parameters fli 
shows the relative contribution of each type of matter to 
the critical energy density at any time. The allowed range 
for any density parameter is then < fii < 1, and fi^ = 1 
means that the evolution of the universe is dominated by 
the i-th fluid. Whenever this happens we usually speak 
of the i-th fluid dominated era in the evolution of the 
universe. 

Also, ^ifi would represent the current relative contri- 
bution of pi^o to the current critical energy density pc,o- 
Hereafter, a subscript '0' will denote current values for 
any quantity. 



E. The curvature of the universe 

There is a normalization issue about the curvature of 
the universe we should be careful with. It is often said 
that the curvature constant can be normalized to take 
three different values: a) k = —1, for an open (Hyper- 
bolic) universe; b) = 0, for a flat (Euclidean) universe; 



and c) fc = 1 for a closed (Spherical) uni verse |3q. For 
further simplicity, one usually finds the suggestion of nor- 
malizing the scale factor too, so that its actual value is 
flO = 1. 

On the other hand, the curvature contributes to the 
energy density of the universe, and the product kr^ in 
metric Q should be dimensionless. Thus, we have two 
options. 

• Dimensionless k and r, and a scale factor a(t) with 
dimensions of length. In this case, we can normal- 
ize the curvature constant as mentioned above, but 
we should notice then that the scale factor can- 
not be normalized arbitrarily as its actual value 
is set by the curvature density parameter as oq = 

• Dimensionless scale factor a(t), r with dimensions 
of length, and k with dimensions of length"^ . Thus, 
the curvature constant is given by fc = —a^HQUkfi- 
Without loss of generality, we can choose to nor- 
malize the scale factor as ao = 1. 

For convenience, the second option will be used through- 
out this paper. 

The constant of curvature is given, at any time, by 
k = —a?VlkH^ ■ This is an interesting relation, since dn = 
H^^, which is called the Hubble length, provides us of 
an estimate of the size of the observable universe at any 
time H nil Hence, k tells us of the deviation 
of the spatial part of metric Q from the flat case. For 
scales r for which y^jA^jr = a^/lTlklH r ^ 1, the universe 
can be considered to have an Euclidean space geometry. 

The current value of the Hubble parameter is Hq = 
TOkms-iMpc-^dimillllll, which implies that the 
current Hubble distance (an estimation of the current 
size of our observable universe) is do, if — 4, 300 Mpc. 
Also, the current value of curvature's density parameter 
is flo,k = 0.0l|2,|^, which seems to indicate we live in 
an (slightly) open universe. 

Thus, the universe we actually see can be plainly taken 
as Euclidean, since the scales at which the (current) cur- 
vature of the universe could be appreciable in the metric 
at the present time (> dnfi/ \/\^kM\ — 43, 000 Mpc) are 
larger than the distance to the farthest object we can 
observe! [33 



III. THE EFFECTIVE GRAVITATIONAL 
POTENTIAL OF THE UNIVERSE 

We will now present the first method to describe the 
expansion of the universe according to the type of matter 
is made of. First, we notice that each energy density can 
be given in terms of the actual value of its corresponding 
density parameter as pi = fli^pco a~^'^% see Eqs. 
and 
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Second, the Friedmann equation Q can be rewritten 
the form|l|l[ll[3 



1 

y + V{a) = -nkfl 



(9) 



where now a dot means derivative with respect to the 
dimensionless time t = H^t. The cosmic time is normal- 
ized in terms of the actual value of the so-called Hubble 
time Hq^ = 14 Gy. The effective gravitational potential 
V{a) explicitly reads 



1 " 



-37,+2 



(10) 



Eq. (|5J resembles the conservation of energy for a par- 
ticle with "space" coordinate a{t) and constant energy 
(l/2)rifc_o|!3|- We see that a flat universe has zero total 
energy, and a closed (open) universe has negative (posi- 
tive) total energy. Also, we would like to stress out that 
the actual values of the density parameters are not all 
independent, but are related through the Friedmann con- 
straint at the present time, 



is only the n-th fluid which has 7„ < 2/3. Using the 
condition V'{ac) = in Eq. (|12a|l . we can write 



^"(«-) = -o E(37. - 2)(7. - ln)^l^,oa-'■'^ . (13) 



All the terms inside the sum are positive by assumption; 
therefore, V"{ac) < 0. The final result is the same if 
more than one fluid have an equation of state 7 < 2/3. 
From this we conclude that any critical point corresponds 
to a maximum. 



IV. DYNAMICAL COSMOLOGICAL SYSTEM 

In the presentation of the second method, we find con- 
venient to take the density parameters as the dynamical 
variables themselves. It is then our purpose in this sec- 
tion to show how to write the Einstein equations as a dy- 
namical autonomous s yste m (for a comprehensive read- 
ing of such systems seep^: and|2ll2ll2ll2ll27ll2^ for 
applications in Cosmology). 



'i,Q 



k,0 ■ 



(11) 



A. The general case 



The gravitational potential H1U|) is negative definite, 
V{a) < 0, if all of the (actual) density parameters are 
positive definite, fli^o ^ 0- This fact may imply in some 
cases the presence of at least one maximum in the po- 
tential. This can be verified by direct calculation of the 
first and second derivatives, 

1 " 

V'{a) = -^(37,-2)a,,oa-3'^'+S (12a) 

i 

1 3 " 

V"ia) = -r(a)-^ V(37.-2)7,f7,.oa-3^-(12b) 
a 2 ^ — ' 

i 

The existence of a critical point is not a trivial thing, 
as it marks the value of the scale factor at which d = 0, 
since Eq. (|10|l is equivalent to the acceleration equation 
d — —V'{a). At a = Gc, the acceleration of the uni- 
verse's expansion vanishes. Therefore, the existence of a 
critical point tells us that the universe should have had a 
decelerated and an accelerated expansion at some stages. 

When does a critical point exist? From Eq. H12a|) , there 
cannot be a critical point if all of the equations of state 
7i > 2/3; the universe will always decelerate in such case. 
In other words, it is necessary the presence of at least one 
perfect fluid with an equation of state with a value less 
than 2/3 for the universe to have an accelerated expan- 
sion. 

Actually, it can be proved in the general case that if 
there is a critical point Oc such that V'{ac) = 0, then 
V"{ac) < 0. For this, let us assume that {n — 1) matter 
fluids have an equation of state 7^ > 2/3, and that it 



Let us assume that there are n+1 perfect fluids present 
in the universe. This is the most general case, as we have 
already mentioned that the curvature itself can be seen as 
a special perfect fluid. The Friedmann constraint helps 
us to reduce in one the number of independent variables, 
i.e., it is only necessary to consider n perfect fluids 40]. 

Choosing the {n + I)-th perfect fluid to be absorbed 
by means of Eq. ((T)), the Einstein equations for the rest 
of the n fluids is given in terms of their corresponding 
density parameters as 



(7, - 7„+i) fli - (7j - 7„+i) 



(14) 



where i, j — 1,2, . . . ,n, and 7ri+i is the barotropic equa- 
tion of state of the (absorbed) (n -I- l)-th perfect fluid. 
The prime denotes derivative with respect to the so-called 
number of e-foldings N = hia. To arrive to Eq. ()14|l . we 
have also made use of Eq. (|3a|l in the form 



H 



2 ^ 



(15) 



Exact solution 



System (|14ll may appear redundant and unnecessary 
at first sight. This is because the behavior of the density 
parameters as functions of the number of e-foldings N 
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can be easily found. By definition (see Eqs. (|Hl), the 
exact solutions of Eqs. (|14|l are 



(16) 



where j = 1, 2, . . . , n + 1. It is not difficult to verify that 
Eqs. (fTljl directly follow from Eqs. and (fT?)|l . 

If one only wants to know the evolution of the den- 
sity parameters, then Eq. Hlt)|) suffices to know all of the 
different stages the universe has gone through during its 
evolution. See sections below for some explicit examples. 



B. Fixed points and stability analysis 

However, the interesting thing to note is that Eq. H14(l 
is a dynamical autonomous system of the form x' = f (x) , 
where x = (fii, $72, • ■ ■ , i^n)- In this respect, Eqs. H14|) 
can be seen as a complementary part of Eq. (jSJ, since 
the latter can tell us which attractor properties can be 
found in the dynamical equations and I^J. 

The critical (fixed) points Xc of the dynamical sys- 
tem ifT^ are found by solving the equations f(xc) = 0. 
There are two obvious solutions. 

• Trivial ({n+l)-th perfect fluid dominated) solution, 
for which Xc = 0, ^n+i = 1, and then the expan- 
sion of the universe is driven by the {n + l)-th per- 
fect fluid. Eq. l(T3|) points out that the universe 
has an accelerated expansion, a > 0, (decelerated 
expansion, a < 0) if 7n+i < 2/3 (7^+1 > 2/3). 

There are two important remarks we should be 
aware of at this point. 

— If the (ri + l)-th perfect fluid is the curvature, 
then 7„+i = 7fe = 2/3, and then the universe 
expands at a constant rate, i.e., a = 0. It 
should be noticed, however, that the existence 
of this critical point is forbidden by the Fried- 
mann constraint {Tj) for the case of a closed 
universe (if all of the density parameters are 
positive definite). 

On the other hand, this trivial solution is per- 
mitted for an open universe, and is better 
known as Milne's universe\(i ll2L l29l| . 

— Milne's model should be distinguished from 
the so-called empty static model^^. Eqs. 
tell us that an empty universe (no perfect fluid 
present, zero curvature) is indeed permitted, 
which will remain static, a{t) = const. 



• j-th perfect fluid dominated solution, for which 
rij = 1, and — = 0; that is, the critical 

points are Xc,i = (1,0,..., 0), Xc,2 = (0,1,..., 0), 
etc. Eq. ifT^ again indicates that the universe will 
accelerate (decelerate) its expansion if the domi- 
nant equation of state is such that jj < 2/3 (7^- > 
2/3), which is the same conclusion we reached at 
in the previous section. 



Once found, the stability of the critical points can be 
established by a first order perturbation analysis[2^, in 
which one considers a small perturbation u in the form 
X = Xc + u. Hence, Eqs. (|14|) are linearized in the form 
u' = Mu, where 



M 



ji 



dfj_ 
dxi 



(17) 



are the elements of the perturbation matrix Ai . 

If the eigenvalues of the matrix M have all negative 
(positive) real parts, then the critical point is called stable 
(unstable). If neither, it is then called a saddle point. 

Back to our case, we have to solve the set of equa- 
tions (|14|1 for Vl'j = 0, and to study the eigenvalues of the 
corresponding perturbation matrix. Let us suppose that 
we want to investigate the stability of the fixed point cor- 
responding to the domination of the i— th perfect fluid, 
Xc,. = (0,0,..., a, = 1,...,0). 

After careful calculations, the elements of the pertur- 
bation matrix are explicitly given by 



Mji = 3(5j7 (7» - 7j) + 3(5j. (7/ - 7n+i) 



(18) 



where 5ji is the Kronecker delta, 5ji — \ li j — I, and 
5ji ~ otherwise. 

The matrix M. is almost a diagonal matrix, except for 
the non-zero elements in the row j = i; however, the 
calculation of its eigenvalues is a simple matter. It can 
be proved that the eigenvalues u>i are solutions to the 
algebraic equation 



nr^V [3(7.-70-^/1 = 0, 



(19) 



where I runs through the permitted (n+l) values, except 
for the the particular value I — i, and then there are only 
n eigenvalues. 

In principle, the trivial critical point should be treated 
separately, and the eigenvalue equation in this case is 



nr=i [3 (7«+i - 7/) - 



0. 



(20) 



Naively, Eq. H20|) seems to be a particular case of Eq. H19|l , 
one in which i' — n + l. However, this is not formally so, 
as the n+ 1-th perfect fluid does not appear in Eqs. H14|l . 

Therefore, we come to the following general conclu- 
sions. 

• The critical point corresponding to the domination 
of the perfect fluid with the largest equation of state 
is unstable. 

• The critical point corresponding to the domination 
of the perfect fluid with the smallest equation of 
state is stable. 

• All other critical points are saddle points. 

The above statements directly follow from Eqs. I|19|) 
and (|20|l . as the positivity or negativity of the pertur- 
bation eigenvalues for a particular case depend on the 
relative value of the dominant equation of state with re- 
spect to the others. 
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V. HOW DOES THE UNIVERSE EXPAND? 

In this section, we will work on particular and simple 
models of the universe, and draw their dynamics accord- 
ing to their material contents and the methods discussed 
in the previous section. 



A. Our universe, or the Concordance Model 

We now turn our attention to the so-called concordance 
model (CM), also known as the ACDM (Lambda Cold 
Dark Matter) model, the model of the universe the cos- 
mological observations altogether seem to favor [lLl^l2ll|. 
It contains radiation (relativistic particles), matter and 
vacuum energies in the proportions fir ft = 10~^, il,n,o = 
0.266, i^Afi = 0.732, and Qk,o = 0.01, respectively. 

For numerical purposes, we take ilr,Q = 10~^, fl,n,() = 
0.3, r^A^o = 0.69, and flk.o = 0.01, so that the Friedmann 
constraint is accurately accomplished at the present time, 
see Eq. ((TTJ- 



1. Effective gravitational potential 
The effective gravitational potential of the CM is 

Via) = -U^ + ^ + n^^oa') . (21) 

As the density parameters are positive definite, we see 
that the gravitational potential has one critical point. 
As radiation does not contribute significantly, the criti- 
cal point is approximately at ac — \/f^m,o/(2^lA.o) = 0.6. 
As said before, this critical point is a maximum^ as 
V"{ac) = — rir^oaj"* — 317a, < 0. For completeness, we 
plot V{a) in Fig.^for a CM universe. 

The curvature of the universe is so small that one 
hardly distinguishes our universe from the 'Flat' case in 
Fig. ^ and the maximum of the effective gravitational 
potential is negative, V(ac) < 0. These two facts imply 
that the CM universe will never stop expanding, had a 
past decelerating stage, and is currently accelerating its 
expansion. 

Notice that, keeping the radiation and matter present 
contributions fixed, the larger the vacuum contribution 
is, the sooner the expansion starts accelerating. 



2. Dynamical system 
Explicitly, the dynamical system of the CM is 

n'^ = Qr (2r2r + - - 2) , (22a) 

= Qrn (2r2r + - '2flA ~ 1) , (22b) 

n'^ = Qa (2f7r + - 2f7A -I- 2) , (22c) 
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FIG. 1; (Top) The effective gravitational potential V{a), see 
Eq. I2H . for the so-called concordance model, according to 
the values of the density parameters Qi obtained from re- 
cent cosmological observations. The solid vertical line goes 
through the maximum of the potential, and marks the time 
at which a = 0. (Bottom) The evolution of the energy density 
parameters Qj as functions of the number of e-foldings A'^, see 
Eqs. Notice that there are radiation, matter and A dom- 

inated stages; though the curvature is not zero, a curvature 
dominated era cannot be achieved. The dashed vertical line 
is at the present value of the (normalized) scale factor ao — 1 
(N = 0) in both figures. 



and its critical points, together with their respective sta- 
bility analysis, are shown in Tabled the eigenvalues of 
the perturbation matrix (|17|l are given in Table UTI 



TABLE L Critical points for the CM dynamical system (1^ . 
As the curvature is positive, the trivial solution does not exist 
for the CM model. 



Domination 




Xo 


Stability 


Radiation 





(1,0,0) 


Unstable 


Matter 





(0,1,0) 


Saddle 


A 





(0,0,1) 


Stable 



As explained before, the trivial solution is also a so- 
lution of the dynamical system (|22|l . but its existence is 
forbidden on physical grounds because of the Friedmann 
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TABLE II: Eigenvalues of the perturbation matrix 1171 cor- 
responding to the critical points shown in Table ^ see also 
Eqs. (EHJ and lO- 



Domination 






UJ3 


Radiation 


1 


2 


4 


Matter 


-1 


1 


3 


A 


-2 


-3 


-4 



constraint. 

The equations of state are 7^ > 7m > 7a, and thus we 
conclude that radiation domination is an unstable point, 
matter domination is a saddle point, and A domination 
is the only stable point. The universe will always reach 
a A dominated era at some epoch, and will remain in 
it thereafter. We want to point out that this is true 
whatever the actual contributions of each component arc. 

The evolution of the density parameters Qj is also 
shown in Fig^ As discussed before, one can see the radi- 
ation, the matter and the A dominated eras, but the cur- 
vature dominated solution is never achieved, even though 
its contribution is noticeable at the transition between 
the matter and A dominated eras. 
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FIG. 2: The effective gravitational potential V{a), see 
Eq. H21^ . for an Einstein's static universe with Qm,o ~ 0.3, 
Qa,o = 1.713, and ^lk,o ~ —1.013; the critical point is at 
ac — 0.44. The above values are determined from the con- 
ditions V{ac) = {l/2)Qk,o, V'{ac) — 0, and the Friedmann 
constraint lllll . 

universe can expand even in the presence of a finely-tuned 
cosmological constant. 



B. Einstein's static universe 

Originally, Einstein considered the universe as static, 
and he introduced a cosmological constant to allow his 
equations Q to have a static solution. That is, Einstein 
found a solution with a = d = 0. For an interesting 
discussion on Einstein's cosmological model see Chs. 2 
and 4 in[^. 

1. Effective gravitational potential 

According to Eqs. © and (|12a|l . Einstein's static so- 
lution corresponds to a universe located at exactly the 
critical point V'{ac) = 0, for which V{ac) = {l/2)i}k,o- 
That is, the universe has just the enough total energy 
to be at the critical point of its effective gravitational 
potential. 

Unfortunately, we have already learned that the crit- 
ical point of V{a) that Einstein found must be a max- 
imum, and therefore is an unstable point. Hence, the 
universe should expand or collapse, but cannot remain 
static! An example of a static universe with density pa- 
rameters ^m,o = 0.3, r^A.o — 1.713, Qk,o = —1.013 is 
shown in Fig. |2| 

Some time after Einstein proposed its static universe, 
the expansion of the universe was discovered, and Ein- 
stein thought he had made a mistake in introducing a cos- 
mological constant to have a static universe; then called 
the latter his 'biggest blunder'. 

It took a bit longer for cosmologists to recognize that 
Einstein's original model is unstable)^ l30l|. and that the 



2. Dynamical system 

The dynamical system analysis proposed in a previ- 
ous section can, in principle, be applied to the Einstein's 
static universe. The conclusion is that the critical points 
are the same as those of the concordance model, see Ta- 
ble and |n] This is because the existence and nature of 
the critical points only depend on the equations of state 
of the material content of the universe. 

However, Einstein's universe is a case we have to deal 
with carefully, because one cannot define a critical den- 
sity at the maximum of the potential V{ac), since the 
latter this time implies the vanishing of the Hubble pa- 
rameter H. 



C. Dark energy 

Since the discovery of the accelerated expansion of the 
universe in the observations of Type la supernovaeQ, 
cosmologists have been wondering whether the energy 
responsible for the acceleration is a cosmological constant 
(7 — 0, p\ — const), or another exotic kind of matter, 
usually dubbed dark energy, with an equation of state 
that lies in the range 2/3 > 7x > and px 7^ const (we 
shall use an X to denote dark energy quantities) 

For concreteness, we will next work on the case of a 
CM universe in which vacuum energy is changed by a 
dark energy component with a constant equation of state 
7a: = 1/3. It should be warned that this is not the most 
general case of dark energy, but we consider a constant 
7a: for pedagogical reasons only. 
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FIG. 3: (Top) The effective gravitational potential V{a) for a 
CM universe in which the vacuum energy is replaced by a dark 
energy component with equation of state 7x = 1/3. (Bottom) 
The corresponding evolution of the density parameters Qj , see 
Eqs. (I16L as functions of the number of e-foldings A''. As in 
Fig.Q there are radiation, matter and dark energy dominated 
stages, without a curvature one. The dashed vertical line is 
at the present value of the (normalized) scale factor ao — 1 
(N = 0) in both figures. 



1. Effective gravitational potential 



The effective gravitational potential and the evolution 
of the density parameters are shown in Fig. |21 As in the 
standard CM model, there is a maximum in the effective 
gravitational potential, but the behavior at late times 
differs from that of the CM. This time the maximum of 
the potential is at Oc — \/^m.,o/^x,o — 0.65. That is, 
it takes longer for this universe to have a dark energy 
dominated, and then accelerated, stage than in the case 
of a cosmological constant. This was to be expected, as 
we mentioned before that a cosmological constant is an 
extreme case for the equation of state of a perfect fluid. 



TABLE III: Critical points for the dynamical system 12311 of 
a universe containing a dark energy component. 



Domination 




Xo 


Stability 


Radiation 





(1,0,0) 


Unstable 


Matter 





(0,1,0) 


Saddle 


X 





(0,0,1) 


Stable 



2. Dynamical system 

The dynamical system for a universe with a dark en- 
ergy component is 

r?; = Qr {2nr + -nx-2] , (23a) 

n'„^ = f7„, [2Qr + f^m -nx-l] , (23b) 

n'x = [2nr + n„, - {fix - 1)] , (23c) 

and the existence and nature of its critical points are 
shown in Table ITTIl 

There is only one stable solution, which is the dark 
energy dominated one. This will make the universe have 
an accelerated expansion at late times. We see that the 
model in this section is similar to the CM, but an accu- 
rate enough measurement could in principle distinguish 
between the two. So far, the CM seems to be the chosen 
one, seeP,|3|. 



TABLE IV: Eigenvalues of the perturbation matrix I17II cor- 
responding to the critical points shown in Table ITTTl 



Domination 




LiJ2 




Radiation 


1 


4 


3 


Matter 


-1 


3 


2 


X 


-3 


-2 


-1 



VI. FINAL COMMENTS 

We stated at the beginning of this manuscript our in- 
tention to understand the solutions of Einstein's equa- 
tions in the case of a homogeneous and isotropic universe. 
After all of the mathematical work done in the previous 
sections, we can extract some general conclusions on how 
the expansion of the universe depend on its material con- 
tent and curvature. 

Assuming all of the energy densities in the universe are 
positive definite, the effective gravitational potential H10() 
is negative definite. Therefore, 

1. An open universe will expand for ever, irrespective 
of the material content and of how large is its neg- 
ative curvature. If there is a component with an 
equation of state such that 7 < 2/3, an open uni- 
verse will enter an accelerated stage at late times. 
Otherwise, it will decelerate but the expansion rate 
will never vanish. 
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2. A flat universe will always expand, and will accel- 
erate its expansion if there is a component with an 
equation of state such that 7 < 2/3. Otherwise, the 
expansion rate will asymptotically vanish (d — *■ as 
t —^ 00). 

3. A closed universe will have an accelerated expan- 
sion if there is a component with an equation of 
state such that 7 < 2/3, and if its curvature is such 
that ^lk,o < 2V{ac). If at least one of the previous 
conditions is not accomplished, then the universe 
will recoUapse. 

Some other comments are in turn. It was shown how 
the dynamics of the universe can be studied using two 
complementary approaches. For instance, the gravita- 
tional potential method gave a clear proof of the insta- 
bility of the Einstein's static universe, whereas the dy- 
namical system method revealed the stability nature of 
the different domination stages our universe has passed 
through. It is our hope that the two methods revised in 
this paper can help undergraduate students to deal with 
the expansion of a homogeneous and isotropic universe. 

Moreover, we would like to stress the fact that such 



methods are taken seriously for the analysis of a great 
variety of cosmological models. Already mentioned was 
the use of the 'gravitational potential' method inj^ 0, 
l3ll |. On the other hand, the use of dynamical systems 
is widely known in the specialized literature; for a non- 
exhaustive hst of examples, see Refs.jHil IH IH HI 
l2^ and references there in. 

In our opinion, it is important to teach undergraduate 
students different methods and techniques currently used 
in specialized research. For we cannot say to what exten- 
sion the methods presented in this work can be used and 
generalized to other equations in Cosmology; but that we 
shall know from the work of future cosmologists. 
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